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Abstract 

We propose a method to measure the neutrino mass kinematicahy using beams of ions which 
undergo beta decay. The idea is to tune the ion beam momentum so that in most decays, the 
electron is forward moving with respect to the beam, and only in decays near the endpoint is the 
electron moving backwards. Then, by counting the backward moving electrons one can observe the 
effect of neutrino mass on the beta spectrum close to the endpoint. In order to reach sensitivities for 
mil < 0.2 eV, it is necessary to control the ion momentum with a precision better than 5p/p < 10~^, 
identify suitable nuclei with low Q- values (in the few to ten keV range), and one must be able to 
observe at least 0(10^^) decays. 



1 



I. INTRODUCTION 



There is now convincing evidence that neutrinos have mass owing to the confirmation of 
neutrino oscillations. Data from atmospheric jl|, ^ , solar [sl, 0, 0] , reactor [g], Q, S] and long- 
baseline experiments j^, [lo| are well described by three- flavor neutrino oscillations, with the 
following values for the oscillation parameters [11|: 
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where Am^j = rrP- — mf are neutrino mass-squared differences, and 6ij are the lepton mixing 
angles. The third mixing angle is only bounded from above with current data: 



sin^^i, < 0.040 



(2a) 
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Neutrino oscillation experiments have yet to determine the sign of Amg^ and so we are left 
with two possible neutrino mass orderings. The case Am^i > is referred to as 'normal 
mass ordering' while the case Amg^^ < is labelled 'inverted mass ordering'. In addition, 
the absolute mass scale is presently unknown. If the absolute scale is small there can be a 
pronounced hierarchy between the masses dependent on the ordering. In particular we can 
have 

1. Normal mass hierarchy (NH): mi <C m2 <C ms 



2. Inverted mass hierarchy (IH): 



m3 <^ mi < m2 



If the absolute scale is large compared to the mass squared splittings, the spectrum is 
referred to as 'quasi- degenerate' (QD), where mi ~ m2 — m^ = m^ with m^ ^ lAmg^l. If 
the lightest neutrino mass is not negligible and of the order y^jAmfTT) then there exists a 
partial hierarchy. 

At present there are three approaches to determining the neutrino mass scale: through 
the measurement of the electron energies near the endpoint of a beta-decay spectrum, mea- 
surement of the half-life in a neutrino-less double beta-decay (if neutrinos are Majorana 
particles), and via study of cosmological data; see Ref. 1^ for a review. An incomplete 
list of phenomenological studies of observables sensitive to absolute neutrino masses can be 



found in Refs. [l3|, [ij, uM, UM, llTI, ll8|, ll9|]. Both the observation of neutrino-less double 
beta decay and cosmology are model dependent approaches. In particular, it is not known 
whether neutrino-less double beta decay exists, and if were observed, there could be other 
mechanisms such as a SU{2)w triplet Higgs [20|, Leptoquarks [2]|, and R-parity violation 
in Supersymmetry [2^ . 

The present bound on the neutrino mass from kinematic studies of beta decay endpoints 



was obtained from the Mainz and Troitsk Tritium beta decay experiments |23l . |2J, |25| 



idpo 
3: 



mf < 2.3 eV 
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where m'^ is the effective neutrino mass in beta decay (see below). The next iteration of 
the Tritium beta decay technologies is the KATRIN experiment, which should be able to 
place a limit m"^ < 0.2 eV [26|. However it is unlikely that these approaches will scale to 
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lower neutrino masses. They can only make a positive identification of the neutrino mass if 
the hierarchy is quasi-degenerate. Rhenium calorimeters, such as proposed by the MARE 



collaboration 27|], are also expected to achieve a limit m'^ < 0.2 eV. 

Thus, new experimental technologies are required to reach the level of hierarchical masses. 
Furthermore, it will be important to verify the mass-mixing hypothesis for neutrino flavor 



conversion, as there are non-mass proposals which can cause neutrino mixing [28|, |29|. This 
means measuring the kinematic effect of not only the absolute mass but the mass differences 
as well. This requires precision at the -y/ Am|-^ ~ 0.05 eV, or even at -\/Am|^ ~ 0.009 eV 
level. One new approach towards the neutrino mass, using ultra-cold atoms has been pro- 



posed recently in [30 



Here we present an idea to search for beta-decay spectrum endpoint distortions using 
radioactive ion beams. By tuning the boost factor of the ions, only electrons very close to 
the endpoint of the beta spectrum move in the direction opposite to the beam direction in 
the laboratory frame. In principle, this allows one to search for kinematic effects of a non- 
zero neutrino mass by counting the electrons with backward trajectories. We explore this 
possibility by performing preliminary sensitivity estimates. We specify the most important 
requirements on the setup in order to achieve sensitivities below 0.2 eV. 

Radioactive ion beams are currently being considered as a possible source of neutrinos for 



a future long baseline neutrino oscillation experiment [3l|. These "beta beams" are subject 



of intense phenomenological and R&D studies. Beta beams with small ion boost factors, 
7 ~ 10, have been discussed in the neutrino literature for cross section and nuclear physics 
measurements at low energies, whilst boosts for long baseline experiments typically range 
7 = 80 — 650. See [3^, l33| for reviews and references. The setup proposed here has many 
aspects in common with such beta beams, though the boost factors needed in our case are 
very close to one: 7 ~ 1 + Q/rrie or v/c ^/2Qjrri^, where the Q- value of the decay is 
assumed to be small compared to the electron mass: Q <^ rrie. In this work we do not 
propose a specific experimental scheme for the a measurement of neutrino mass. Whether 
such an experiment can be integrated in a "high-7" beta beam facility is an interesting 
question to be addressed in future studies. 

This proposal has significant challenges that must be considered. These include the 
identification of ions with low Q- values in order to maximize the effect of the neutrino mass 
{Q in the range of few to ten keV), as well as the required number of useful decays of order 
^ 10^^. Furthermore, the momentum spread of the ions in the beam has to be less than 
10~^. This can be achieved either with classical cooling techniques or by exploring the use 



of "crystallized beam" technology [3J, |35 



This paper is structured as follows. In section [TTl we briefly review the phenomenology of 
the measurement of the neutrino mass using endpoint studies. In section UTTl we outline our 
experimental proposal, while in section HVl we present simulations of the precision this setup 
can achieve, and the requirements on various ingredients to reach m'^ < 0.2 eV. In section IVl 
we comment on challenges of the proposed measurement. A summary follows in section IVll 
In Appendix |X] we provide supplementary information on our measurement strategy. 
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II. BETA DECAY ENDPOINT PHENOMENOLOGY 



The standard approach for a direct mass measurement is the analysis of the endpoint 



region of a beta decay 23|, |2J, |25j, see Ref. 36| for a review. A measurement is made by 



reconstructing of the electron spectrum of the decays of ion /, 

I — ^I' + e- + Ue- (2.1) 

It should be noted that such a measurement is sensitive to the mass scale of the anti-neutrino 
which is assumed to be identical to its neutrino counterpart owing to CPT invariance, but 



it need not be [28|. Nuclei which decay through positron emission also have competing 
electron capture decay channels. Electron capture dominates for proton-rich nuclei with low 
Q-values. Positron decay is kinematically forbidden for Q < 2me effectively prohibiting the 
measurement of the neutrino mass by this method. 

The three mixing angles and CP-phase parameterize a mixing matrix, with elements 
Uaii which relate the neutrino mass eigenstates with the eigenstates that participate in the 
weak interaction. This mixing of the neutrino mass eigenstates means that the spectrum 
of the electrons from the decay should be considered as the incoherent sum of the spectra 
associated with each neutrino mass eigenstate, 

where 



PpEp{E^^^ - Efs)JiE^.,, - Epf - mf F{Z, Ep) S{Ep) [1 + 5r{Z, Er)] (2.3) 



dEp 

are the individual electron spectra and mi is the mass for eigenstate i. A recent evaluation 



of the spectrum can be found in |37l . l38l |. Here Ep is the kinetic energy of the electron and 
-Emax is the maximum electron kinetic energy for zero neutrino mass = ("endpoint 
energy"). E^ax is given by 

^max = -rrie^Q, (2.4) 

where M and M' are the masses of the mother and daughter ions / and /', respectively, and 
we define the Q- value to be Q = M — M' — mg.^ The approximate relation -Emax ~ Q holds 
under the assumptions (M — M')/M <^ 1 and mg <^ M. Further, S{Ep) is a form factor 
that contains the nuclear matrix element and constants 

S{Ep) = Gl cos^e \Mp{Ep)\\ (2.5) 



Strictly speaking the Q-value is defined as the mass difference of neutral atoms. Here we denote with M 
and M' the actual masses of the ion before and after the beta decay, and therefore we explicitly include 
the electron mass in the expression for Q. This holds up to corrections of order of the binding energy of 
the electron. 
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The Fermi function F(Z, Ep) describes the Coulomb interactions of the ejected particle. For 
nuclear radius R = 1.2A^^^ and the definitions r] = aZEp/pp and 7 = (1 — (aZ)^)^/^, it 
takes the form 

V ft / '|r(27+l)l 1 - exp(-27rr,) ^ 

The Coulomb correction increases the decay rate for (3~ decays and decreases it for (5^ decays 
because the attraction to the nucleus enhances the decays with low momentum. 5r is the 



contribution from electromagnetic radiative corrections [37j which are usually negligible. 



The electron spectrum in Eq. 12.21 is parameterized by the 3 neutrino masses and the 
mixing angles ^12,^13- However, for experiments with resolution worse than \/| A mjil, we 
may parameterize the spectrum with a single effective mass Trf^ , see e.g., 39|, |40|, |41|: 



m 



(normal ordering) 
V "^min + I ^"^31 1 +(9( v^Am2i) (inverted ordering) 



(2.7) 

where m^^^ is the mass of the lightest neutrino, and the approximate expressions hold up 
to terms of order a/ ISm^^ ~ sin ^13 A/]Am|^ ~ 0.01 eV. Then the spectrum becomes 



dE, 



oc ppEp{E^^^ - Efs)J{E^,^ - Epf - (mf )2 F(Z, Ep) . {21 



In Fig.dl we show the effective mass Tjf^ as a function of the lightest neutrino mass mmin- 
The behavior for normal mass ordering is shown in red whilst the blue lines are for inverted 
mass ordering. For a small minimum neutrino mass, the eigenstates separated by the solar 
splitting determine the size of the effective mass. For normal ordering, this pair is positioned 
at mmin resulting in a small effective mass. For inverted mass ordering, however, the solar 
pair is separated from mmin by the atmospheric mass splitting. The effective mass is a factor 
5 larger as a consequence. For mmin ~ a/ Am^^^ = 0.049 eV, the atmospheric splitting is 
not dominant and the effective mass does not discriminate between the two orderings. 

The goal of future absolute neutrino mass experiments is to push the sensitivity below 
0.04 eV. Failure to measure the neutrino mass above this level identifies the mass hierarchy 



to be normal, on the assumption of neutrino mixing [42]. Failure to measure neutrino mass 
above rrfj^ ~ 0.006 eV would be inconsistent with our current understanding of neutrino 
mass and mixing given in Eq. 11.11 While reaching such sensitivities would be the ultimate 
goal of neutrino mass measurements, in this work we are slightly less ambitious. We have 
in mind experiments with sensitivities in the region 0.04 eV ^ TrfJ^ ^ 0.2 eV, and hence we 
can describe the spectrum by Eq. 12.81 using the effective neutrino mass m^J^ . 

From the expression for the spectrum in Eq. 12.81 it becomes clear that the effect of 
the neutrino mass is larger for decays with a small Q-value. Most previous and present 
experiments are using Tritium, with a relatively low endpoint energy of 18.6 keV and a half 
life of 12.3 y. Alternatively, the MARE project uses ^^^Re, with an even smaller endpoint of 
2.47 keV, at the price of a much longer half life of 43.2 Gy. In this work we do not choose a 
specific ion. Instead we perform the analysis as a function of the Q-value. The identification 
of a suitable ion is central if the approach proposed in this paper is to be realized. 
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FIG. 1: The effective mass as a function of lightest neutrino mass TOmin for normal mass ordering (red) and 
inverted mass ordering (blue). In all cases, ^13 = 0. The solid lines have been simulated for the 3cr upper 
bound for the solar and atmospheric parameters as defined in whilst the dashed lines use the 3cr lower 
bounds. 



III. THE CONCEPT 



In this article we consider the possibility of using a very low boost ion beam as a tool to 
measure the neutrino mass scale, by observing backward moving electrons in the laboratory 
frame. A low boost [y/c ~ ^/iQjm^ radioactive ion beam is sent through an evacuated 
chamber with a weak magnetic field parallel to the beam line. A detector is set up on the 
back wall of the chamber to record the number of electrons still traveling backwards after 
the boost, see Fig. [21 Therefore, the purpose of the boost is to perform a cut on the electron 
momenta, only selecting electrons very close to the spectrum endpoint. 

Let a boosted ion have a velocity vi in the laboratory frame, and an electron have 
velocity f e in the rest frame of the ion. Further, let Q be the angle between v^. and the beam 
direction. Hence, f I = cos Q is the electron velocity component parallel to the ion beam. 
Any electrons that satisfy + < will appear in the backward direction in the laboratory 
frame. Let the maximum momentum of the electrons in the ion rest frame be Pmax? with 

Pmax = V^max(^max + 2me) ^ ^j2m^Q, (3.1) 

where the approximation holds for small Q- values: -Emax ~ Q <C mg. Note that in this 
regime, ions, as well as electrons, are non-relativistic and we can use Galilean velocity 
addition.^ Let us define an ion momentum p^? such that an electron with momentum Pmax 



^ We have checked explicitly that rclativistic effects are small and can be neglected to good approximation. 
Note that the neutrino is rclativistic and its momentum is not Galilean invariant. 
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electron detector 




ion beam 



FIG. 2: Diagram of the proposed experiment. The ion beam enters an evacuated cavity whose back wall 
holds an electron detector. Each ejected electron follows a helical trajectory. Electrons moving in the 
backward direction in the laboratory frame are counted by the detector. Sec also Fig. [71 



emitted in the backward direction is at rest: 



M 



nip. 



and 




0.2 M 



Q 



lOkeV 



1/2 



(3.2) 



From this relation we obtain the boost factor of the ions in the beam. We have 7 = E'j/M 
^/l + 2Q/me, which is close to one for the case of interest, when Q -C nie. Now one can 
alter the ion boost slightly from that value by a small Ap/ = p^j — pi such that only electrons 
close to the endpoint will go in the backward direction. More precisely, the ion boost will 
perform a cut such that only electrons with parallel momentum satisfying 



Pmax ^ ^ P\\ ^ Pmax 

will travel in the backward direction in the laboratory frame, where from Eq. 



(3.3) 
we get 
(3.4) 



Hence one selects electrons with momenta close to the endpoint. Note that fromEq. 13.41 
the cut in electron momentum e is related to the change in the ion momentum by the tiny 
number irie/M. Combining this relation with Eq. 13.11 one can translate e into a cut on 
electron energy: 




Api ^ 2Q 



Api 

P'l 



(3.5) 
0.1 eV 



me M 

Therefore, assuming Q — 5 keV, sensitivity to neutrino masses of order Ai?^ 
requires control of the ion momentum at the level of 5p/p ~ 10~^. 

The total number of electrons going in the backward direction can be calculated by 
transforming to cylindrical coordinates: p^dp_[_dp\\ = 2p{Ep + me)dEp with p^ = p]_+p'\\ = 
Ep{Ep + 2me). Then one finds 
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FIG. 3: The fraction of electrons that travel in the backward direction (Eq l3.6p in the laboratory frame as 
a function of the cut on the electron momentum parallel to the beam, e, sec Eq. 13.31 A Q-value of 5 keV 
has been considered for a selection of neutrino masses from zero to ~ 1 eV. 

In Fig. [3] we show the number of electrons going backward relative to the total number of 
electrons, which is given by 

iVtot = 2iV(e = p^ax) = / dE^TET- (3-7) 

Jo dEp 

The difference between rrf^ = and m^^ ~ 0.2 eV is at the level of 10"^^ and we conclude 
that in order to reach these sensitivities to the neutrino masses, of order > 10^^ ion decays 
are needed. This number, of course, represents a major challenge for the proposed technique. 

Now we assume a magnetic field B parallel to the beam axis. The electrons will spiral 
with gyro-radius p±/qB {q = 1 for electrons). Hence, electrons with a given orthogonal 
momentum component will hit either the front or back wall of the cavern at a distance from 
the beam axis ranging from zero to -Rmax(p±) = '2p±/B. The distribution of events as a 
function of R encodes some spectral information, the use of which slightly improves our 
sensitivity, as described in Appendix [A] 

A traditional endpoint search, such as KATRIN uses an electric field to select electrons 
close to the endpoint; a cut on the momentum p > pcut such that a shell of the momentum 
sphere in the ion rest frame is selected. In our case we cut on the momentum component 
parallel to the beam. Hence, we use only a slice of the momentum sphere, which has the 
disadvantage that many electrons close to the endpoint are lost. Roughly, the fraction of 
electrons in the shce (p max - e < p\\ < Pmax) and the shell (pcut = Pmax-e < p < Pmi,^) is given 
by (p™'^'^/Pmax)^ ~ e/v^TTieQ ~ 10""^, whcrc p^^"^ ~ V^Phna^ and the value 10~^ is obtained 
for e ~ 5 eV and Q 5 keV. This small factor of useful decays has to be compensated by 
the total number of ion decays, which is part of the reason our proposal requires > 10^^ 
decays. 
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IV. SIMULATIONS AND SENSITIVITY ESTIMATIONS 



A. Description of the analysis 



To extract the neutrino mass from an experiment as suggested here, one would perform a 
measurement at many values of e in order to sample the spectrum shown in Fig. [31 Since the 
count rate increases very fast with e one would optimize the measurement time at a given e 
such that most time is spent close to the endpoint at small e, whereas the necessary time in 
order to accumulate enough events decreases fast with increasing e. For each experimental 
run at fixed e, one can perform a fit to the distribution of events as a function of R, the 
distance of the detected event from the beam axis, as discussed in Appendix \M Before 
commenting on such a fuU-fiedged analysis, we discuss first a simplified analysis using only 
the total number of counts (no radial information) at two values of e. 

At least two data points are needed to simultaneously extract the Q-value of the ion. 
Typical uncertainties on Q are too large compared to the precision required in order to use 
the Q-value as an input for the endpoint measurement [43| . For example, the Q- value for 
Tritium decay has been determined as 18589.8 ±1.2 eV [4^, to be compared to the 0.2 eV 
sensitivity goal of KATRIN. Future high precision measurements at the MPIK/UW-PTMS 
in Heidelberg 45|] aim at a precision for Q of 30 meV, which may be used as cross check 
for the KATRIN experiment. In gen eral it is therefore necessary to fit for the Q- value, in 



addition to the neutrino mass [26|, |43|. Even though the neutrino mass and Q- value are not 
degenerate if full spectral information is available, in an analysis of the total counts at a 
single e, the effect of a non-zero neutrino mass can mimic a slightly larger Q- value. Hence, a 
measurement at a single value for the ion boost is not sufficient to fit for the neutrino mass. 

In the simulations carried out in this study, we combine two experimental runs: one with 
small e close to the endpoint and one with large e. At large e the total count rate is several 
orders of magnitude larger than the change invoked by a non-zero neutrino mass. Hence, 
a 2-parameter fit for the neutrino mass and Q-value at large e will be largely independent 
of the neutrino mass, i.e. one effectively makes a measurement of the Q-value. For small e, 
however, the total count rate becomes comparable to the reduction for non-zero neutrino 
mass. There is a strong neutrino mass dependence in this case which, when combined with 
the Q-value measurement from the run with large e, constrains the neutrino mass. 

Quantitatively, in our numerical analyses we calculate the following 



X = mm 



Our definition is based on Poisson statistics: 



V 



(4.1) 



De, + D,, In 



where 



(4.2) 



iV,, = il + v) Tjmf , AQ) , A. = Tjmr , AQ = 0) . (4.3) 

T^.{m'^,AQ) is the predicted number of events for displacement for a certain neutrino 
mass hypothesis m'^ and a shift in the Q-value of AQ. The pull parameter t] takes into 
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No systematics 



Error on flux: 2 % 




FIG. 4: We show the 90%, 95% and 99% confidence level contours for two degrees of freedom. We assume 
a Q-value of 5 keV, 10^® useful ion decays, and a "true" neutrino mass of = 0.1 eV. Results are shown 
separately for e = 5 eV and e = 100 eV, and their combination. In the left panel no systematical errors are 
assumed, whereas in the right panel, we introduce a 2% error on the flux normalization. No backgrounds 
have been included in these analyses. 

account a systematic uncertainty on the overall normalization of the number of events, 
o"norm, which we assume to be correlated between the two runs. D^^ is the simulated data 
for displacement at an assumed "true value" for m'^ and at the "true" Q-value. The 
final is found by adding the penalty term (///(Tnorm)^ and then minimizing with respect 
to the pull. Owing to the large uncertainties on the Q-value, we always treat AQ as a free 
parameter. 

In Fig. H] we demonstrate this approach for the case Q = 5 keV and m^J^'^ = 0.1 eV, both 
for no systematics and a 2% error on the normalization of the flux. Runs with e = 5 eV 
and 100 eV have been taken with measurement times in the ratio 99:1 such that the total 
number of useful ion decays is 10^^. It is evident from Fig. H] that for large e, the fit only has 
a very slight dependence on the neutrino mass and provides a measurement of Q, where the 
accuracy on the Q-value is sensitive to uncertainty on the normalization of the number of 
useful decays. For small e, on the other hand, we see a strong dependence on the neutrino 
mass which shows little variation with the normalization of the flux. Obviously, a low e 
measurement alone cannot constrain m'^ because of the correlation with AQ. However, the 
combination of small and large e runs provides sensitivity to the neutrino mass. The effect 
of the systematics on the neutrino mass sensitivity is felt through the uncertainty on the 
Q- value. 

We have checked that the two-e run analysis provides already a good sensitivity. Certainly 
there is room for improvement by exploring many runs at different e values with optimized 
measurement times at each position. While a detailed optimization along these lines is 
beyond the scope of the present work, preliminary studies indicate that the sensitivity may 
be improved by a factor of order 10%, similar to exploring the radial information as described 
Appendix [XI 

Fig. m indicates that the proposed measurement may provide a determination of the Q- 



10 



o 

ON 



0.01 
10 





rate only @ £ = 5 and 100 eV 

/^-distribution @ £ = 5, 100 eV 














a 

nonn 



16 



J I 



J I 



J I 



J I I 



10 



17 



J8 



10 

# of total decays 



10 



19 



10 



20 



FIG. 5: Sensitivity to tlic neutrino mass at 90% confidence level as a function of useful decays for Q = 
2, 4, 8 keV. We show the total rate analyses with two ion boosts corresponding e = 5 and 100 eV (solid) and 
an analysis with using also the radial distribution with 20 bins of equal width in R (dashed). 

value at the sub-eV level. Let us stress that this would rely on a perfect knowledge of the 
momentum of the ion in the beam pj. It turns out that AQ is fully correlated with pj and 
what actually is measured is the sum of the two. Therefore, if we are interested mainly in 
a measurement of the neutrino mass, without the ambition to determine also Q, we do not 
need to know pj with very good accuracy, since we simply fit for the sum pi + Q. However, 
it is very important that the momentum spread 6pi of the beam is small enough, such that 
the spectral distortion introduced by m'^ is not washed out. We are going to quantify this 
requirement below. 



B. Requirements to reach sub-0.2 eV sensitivity on the neutrino mass 

In Fig. we present the upper bound on m'^ at 90% confidence level, which can be 
obtained if the true value is m^J'^'^ = 0, as a function of useful decays for Q- values of 2, 4 and 
8 keV. Runs of e = 5 eV and e = 100 eV are considered in the ratio 99:1 such that the total 
useful number of decays sum up to the value shown on the horizontal axis. No backgrounds 
or systematics have been included, and we neglect the momentum spread of the ions. It is 
seen from Fig. \5\ that the approach adopted in this paper places very strong requirements 
on Q and the luminosity. To match the expected sensitivity of the KATRIN and MARE 
experiments of m'^ < 0.2 eV, one requires 4 ■ 10^^, 5 ■ 10^''', 8 ■ 10^^ decays for Q-values of 2, 
4, 8 keV, respectively. A factor 5 improvement down to m'^ < 0.04 eV, which will separate 
the normal and inverted mass hierarchy regions in Fig. [H requires in excess of 10^^ counts 
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FIG. 6: Sensitivity to at 90% confidence level for (5 = 5 keV and 10^^ useful decays (dashed/blue) and 
Q = 3 keV and 10^^ useful decays (solid/red). In the left panel, the effect of the normalization uncertainty 
on the flux is considered; background levels are varied in the center panel; whilst the effect of the momentum 
spread of the initial ion is taken into account in the right panel. 

across the run of the experiment. 

We also show in Fig. [5] an analysis that takes into account the radial distribution of 
the backward moving electrons, as described in Appendix |Al 20 bins of equal radial width 
were used with the same run parameters as the total rates analysis. There is only minor 
improvement when radial information is included in the very high luminosity region, because 
this approach does not reconstruct the electron energy spectrum itself, as electrons with 
transverse momentum p± can strike the detector at any radius less than twice their gyro- 
radius Tg = 2p±/B. The spectral information is thus significantly smeared out and cannot 
be reconstructed accurately. Although neither analysis has been rigorously optimized, one 
does not expect order of magnitude increases in sensitivity. Therefore, in the following, we 
only consider the two-e run analysis described in section IIV A[ 

In Fig. [6] we examine the behavior of neutrino mass sensitivity with systematics, back- 
grounds and the ion beam momentum spread. We consider the cases Q = 3 keV with 10^^ 
useful decays ("high sensitivity") and Q = 5 keV with 10^*^ useful decays ("low sensitivity"), 
and, as before, in both cases we run at e = 5 and 100 eV with a ratio of 99:1. The "low" and 
"high" sensitivity configurations provide nominal sensitivities of 0.21 and 0.071 eV, respec- 
tively. Now we include each one of the three above mentioned effects separately, in order to 
investigate at which level the sensitivity starts to deviate from these idealized numbers. 

From the left panel of Fig. El it is seen that the sensitivity starts to deteriorate rapidly for 
normalization uncertainties greater than 1%. This can be understood in conjunction with 
Fig. m where it was shown how the error on the normalization on the number of useful decays 
affects the overall sensitivity. The systematics severely limits the ability of the e = 100 eV 
run to constrain the Q- value which in turn worsens the sensitivity to the neutrino mass. This 
is true for both cases considered, but especially the high sensitivity setup which becomes 
independent of o"norm only below 0.1%. Naively such an accuracy seems not very difficult to 
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achieve. Apart from traditional beam flux determinations, one could measure the electrons 
emitted in the forward direction by installing an electron detector at the front wall of the 
chamber. The (9(10^^) electrons may provide a means to determine the total number of 
decays with negligible statistical error. 

In the center panel of Fig. [6l the effect of including backgrounds is shown. We assume 
that the magnetic field is adjusted such that the total area illuminated by the electrons 
on the back wall is the same for both the small and large e. A constant background rate 
is taken. The number of background events shown on the horizontal axis is the combined 
background from both runs. Here, the behavior is as one would expect; the low sensitivity 
setup suffers badly with loss in sensitivity even for background rates less than 100 across 
the entirety of the experiment. The high sensitivity setup performs much better; it is able 
to tolerate up to 1000 background events before the sensitivity starts to diminish. 

In the proposed experiment there are several potentially dangerous sources of back- 
grounds. Even a tiny fraction of the 10^*^ to 10^^ electrons supposed to go in the forward 
direction, which could be deflected by some effect to the back wall, would constitute a serious 
background for the measurement. Note that the forward-going decay electrons are high en- 
ergy, and can eject electrons from the forward wall. These could then spiral backwards and 
hit the rear wall. This could be mitigated somewhat by removing the magnetic field in the 
forward direction, and using timing information of the beam bunch relative to the electron 
detection. We also ignore the electric and magnetic field generated by the beam, which will 
affect the trajectory of the ejected electrons. Beam-related backgrounds can be measured 
by adjusting the ion's momentum so that e < and no electrons should have a backward 
trajectory, while the beam-unrelated background can be determined from beam-off periods. 

In the right panel of Fig. El we present the sensitivity as a function of the spread of the 
ion momentum in the beam direction. We assume a Gaussian momentum distribution in 
the beam direction with width 5pi. From Eq. 13. 5^ this translates into a spread on e of 

5e=^2Q^,^-^. (4.4) 
Pi 

In order to include this uncertainty in the analysis we fold the number of events N{e), given 
in Eq. 13.61 with a Gaussian with width 6e. This procedure takes into account the momentum 
spread parallel to the beam, whereas the momentum uncertainty perpendicular to the beam 
is not included. This will be important in reality, especially in a setup that uses the radial 
distribution of electrons. 

In the right panel of Fig. [6] one can see that for the low sensitivity setup, the sensitivity 
to the neutrino mass starts to abate at 6pi/pi ~ 2- 10~^. This is consistent with the relation 
between the uncertainty on the electron energy and the ion momentum spread: 

. (4.5) 

For 6Efs = 0.22 eV and Q = 5 keV, one requires an uncertainty on the momentum to be 
Spi/pi < 2.2 ■ 10~^ for ^ Q. The same estimate for the high sensitivity setup predicts 
6Ei3 ~ 0.075 eV for 6pi/pi ~ 1.2- 10"^. However, from Fig. [6]we observe that the sensitivity 
starts to deteriorate only at 6pj/pi ~ 1.0 ■ 10~^ where, according to Eq. 14. 5^ 6Ep ~ 0.6 eV. 



SE, 



2E,^J^ + O 
Pi 
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This behaviour indicates that the higher count rates can compensate for less precision on the 
ion beam momentum. In summary, to reach interesting sensitivities to the neutrino mass, 
momentum spreads in the range 10~^ to 10~^ are required. 



V. CHALLENGES AND TECHNOLOGY REQUIREMENTS 
A. Total number of decays 

Our analysis indicates that we would require 10^^ — 10^° useful decays to reach sub 



KATRIN sensitivities. Such production is in line with projected EURISOL intensities 46 



47j. However, to store such a large number of ions is a major challenge, especially as they 
will be long lived and possibly have high charge. Space charge effects, the defocussing of 
the beam due to beam-particle interactions, will be significant at non-relativistic velocities. 



B. Momentum spread requirements and coulomb crystals in storage rings 

The mass measurements discussed in this paper can only be realized if the energy spread 
of the ions in the ring is very small. For an electron with kinetic energy Ef^, an uncertainty on 
the ion momentum Spi/pi translates to an uncertainty on the electron energy via Eq. 14. 5[ 
For example, an endpoint electron from an ion with Q = 5 keV, an uncertainty on its 
kinetic energy of 0.1 eV requires Spi/pi = 1 ■ 10^^. The numerical results presented in 
section IIVBI indicate that in order to reach interesting sensitivities to the neutrino mass, 
ion momentum spreads in the range 10~^ to 10~^ are required. Values in this range may be 
achieved with classical beam cooling for low energy ion beams with electron cooling or laser 
cooling. Therefore, if ions with very low Q-values ((9(1 keV)) are available, classical cooling 
is sufficient. Availability of such low Q-values in nature is minimal, however, those that do 
exist are not practical. For example, ^^''Re has Q = 2.6 keV, but it also has a half-life of 
4.5 ■ 10^° years. New coolings techniques will be required for the likely much larger Q-values. 

For cooled low intensity ion beams of 5000 — 10000 ions, a distinct transition to much lower 
momentum spread has been observed with increasing electron cooling current at NAP-M in 



Novosibirsk [48[, ESR at GSI [49| and CRYRING at MSL [50[. This has been interpreted 
as an ordering of the ions in a regime where the energy spread of the beam is too small to 
permit any individual ion to overtake or drop behind its neighbors. Effectively, the ions will 
create a one dimensional string of ions with a minimum distance separating individual ions. 
This distance is defined classically by setting kinetic and potential energy equal, giving 



dmin = {ZqY' / ATTe^kT where kT is the ion temperature 50|. The momentum spread is 
measured to be smaller than 10"^ with the upper limit given by the known ripple of the 
power supplies controlling the electromagnetic fields which confine the beam. The ultimate 
momentum spread is set by the cooling forces acting on the beam. First studies have shown 
that bunching of ordered beam should be possible enabling at least moderate acceleration 



5l[ . To create much higher energy ordered beams it is necessary to develop cooling schemes 
for highly relativistic beams as the ordered beam bunch is unlikely to survive transition. 
Development of electron cooling schemes for ultra relativistic ion beams are underway at 



BNL 5^, for example. The new HE-storage rings at the FAIR facility in Germany will in 
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time enable experiments with laser cooling of highly relativistic beams 53 . 

For the low energy regime, experimental work is in progress in Japan at the S-LSR ring 
in Kyoto 5J]. They have demonstrated one dimensional ordering of protons 55|] and are 



working on special tapered cooling schemes in a dispersion free storage ring for experiments 
aiming to achieve three dimensional ordering of highly charged ions. 

Potential new applications such as the one proposed in this paper could motivate R&D 
on ordered beams. Such efforts could eventually yield higher intensity and/or higher energy 
ordered beams. Studies are needed on increasing the number of ions, and ensuring the 
desired properties of the beam at high densities. All experiments on crystallized beams 



carried out thus far 35 have used less than 10000 ions. 



C. Comments on Q- values 

In this paper, we have established the necessity of > 10^^ decays and Q-values ~ 1 — 10 keV 
if accelerated ions are to be used to measure the neutrino mass at sub-KATRIN levels. This 
measurement therefore needs a low Q- value ion with a short half-life. Rhenium-187 possesses 
the smallest Q-value known: 2.6 KeV; however, the 5/2"*^ — > l/2~ transition has a half-life 
of 4.5 ■ lO^'' y. For the set of beta emitting nuclei with half-lives < 10 y, Ruthenium-106 
and Radium-228 have the smallest Q-values: 39.4 keV and 45.9 keV, respectively. The 
total number of decays will need to be increased by 3-4 orders of magnitude to match the 
sensitivities discussed in Sec. lIVBl if these ions were to be used. 

The most promising ion is therefore Tritium with Q = 18.6 keV, however electron cooling 
is necessary rather than laser cooling because it has no hyperfine structure when ionized. 
Tritium has the added advantage that the number of ions that can be obtained is orders of 
magnitude larger than those discussed in this paper, as it is naturally occurring and does 
not need to be produced in an accelerator. 



VI. SUMMARY AND CONCLUSIONS 

We have discussed the possibility to use radioactive ion beams to study the kinematic 
effects of a non-zero neutrino mass close to the endpoint of the beta decay spectrum. The 
underlying idea is to adjust the boost factor of the ions in such a way that only electrons close 
to the endpoint will have backward trajectories in the laboratory frame. We have discussed 
some requirements of the proposed setup in order to exceed the sensitivity of 0.2 eV of 
the latest generation of Tritium and Rhenium decay experiments. A crucial question is 
whether it will be possible to accelerate enough ions within reasonable time such that of 
order 10^^ — 10^'^ decays can be observed. This issue is also related to the identification of 
a suitable ion with a low enough Q-value (in order to maximize the effect of the neutrino 
mass) , and a small enough half life (in order to have a high enough decay rate) . 

As an example, for a very low Q-value of 2 keV, one needs 4 ■ 10^^ decays to obtain 
the KATRIN sensitivity of 0.2 eV, while 10^^ decays will allow for a m'^ measurement at 
0.04 eV. On the other hand, if no suitable ion with such a low Q-value can be identified 
the requirements on the total number of decays increases drastically: for Q = A{8) keV 
the 0.2 eV sensitivity is reached for 5 ■ 10^^ (8 ■ 10^^) decays, respectively. The sensitivity 
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goal of m'^ < 0.04 eV, which will separate the normal and inverted neutrino mass hierarchy 
regions, requires in excess of 10^^ counts across the run of the experiment together with a 
Q- value of 2 keV. 

Our proposal relies on a very small momentum spread of the ions, at the level of Spi/pi ^ 
10"^. Momentum spreads in the range 10~^ to 10^^ may be reached by classical beam cooling 
techniques such as electron or laser cooling. For even smaller momentum spreads one may 



explore the use of "crystallized beams" [3J, |35|]. This is a phenomenon which has been 



observed for cooled low intensity ion beams, resulting in a transition to an ordered state 
of the beam with momentum spreads of less than 10~^. It has yet to be proven whether 
this technique allows for the extremely high beam intensities required for the neutrino mass 
measurement proposed here. In addition one needs to separate forward and backward moving 
electrons with an precision of order 10~^^. In general forward moving electrons will create 
backward moving electrons by collisions with the walls of the cavern. A suitable experimental 
configuration to overcome this challenge must be identified. 

Given the utmost importance of establishing the absolute value of the neutrino mass in 
a model independent way, we feel that all possible directions have to be explored. We hope 
that our work will stimulate more intense investigations on the possibility to use ion beams 
to pursue this fundamental issue in neutrino physics. 
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APPENDIX A: RADIAL DISTRIBUTION 

In this appendix we comment on the radial distribution of the electron events on the back 
wall of the chamber. One could use a silicon pixel detector which gives position information. 
One can bin the events as a function of the distance from the beam and perform a fit to the 
R distribution for each experimental run at given momentum cut e. 

Consider an electron with charge q = 1 and momentum p_i_ perpendicular to a magnetic 
field B, which is parallel to the beam axis. The trajectory of the particle will be a helix 
with the maximum distance from the beam axis given by twice the gyro-radius 

-Rmax = = — — - . (Al) 
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FIG. 7: Trajectory of an electron in the plane perpendicular to the beam. 



The trajectory of the electron in the plane orthogonal to the beam is shown in Fig. [71 where 
R is the distance from the beam. The electrons will hit the detector with a flat distribution 
in the angle (f. We obtain for the differential count rate (compare Eq. 13.61) : 

From Fig. [7] one finds 

R = 2rg sin — and dip = dR . (A3) 



RB 
2pj 



Hence, we obtain 



dN = — — , dpiidp^ dR . A4) 

Att dEp p{Ep + m,) / _ /^^2 ^ > 

Integration of the momentum for a given e yields the R distribution of the events: 

dN B r'^'" , 1 /•pr''(P^) dV 1 

= — / dpx^ — — / dpw , (A5) 

dR Att ipmin(^) 1^ TrbV . dEfs p{Ep + me) 

2pi 



where p™^^ = ^y2p^y^s,^e — is the maximum perpendicular momentum for given e, p™™(i?) = 
BR/ 2 is the minimum perpendicular momentum required to hit at a distance R from the 
beam, and p^^{p±) = -\/Pmax ~ pI- Note that integrating Eq. IA5I over R from zero to -Rmax 
returns Eq. 13.61 

In Fig. [8] we show some examples for the R distribution. Unfortunately the shape of this 
distribution carries only limited information on the neutrino mass. The reason is that there 
is a significant averaging of the endpoint region. Note that the electrons with p™'^^, which 
have the maximal momentum Pmax, will hit the wall at all radii from zero to -Rmax- Hence, 
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FIG. 8: Number of events as a function of the distance R from the beam axis for 10^* decays at an ion boost 
corresponding to e = 5 eV. The distribution is shown for Q = 3 and 5 keV and for m°*^ = 0, 0.2, 0.4 eV. 

for small R electrons with rather broad range of momenta contribute. Because of this the R 
distribution provides only limited additional information, beyond just the decrease of event 
number with increasing m'^f. This explains the only modest improvement of the sensitivity 
seen in Fig. [51 
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